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Abstract 



We exploit an arbitrary extrinsic time foliation of spacetime to solve the 
constraints in spherically symmetric general relativity. Among such foliations 
there is a one parameter family, linear and homogeneous in the extrinsic 
curvature, which permit the momentum constraint to be solved exactly. This 
family includes, as special cases, the extrinsic time gauges that have been 
exploited in the past. These foliations have the property that the extrinsic 
curvature is spacelike with respect to the the spherically symmetric superspace 
metric. What is remarkable is that the linearity can be relaxed at no essential 
extra cost which permits us to isolate a large non - pathological dense subset 
of all extrinsic time foliations. We identify properties of solutions which are 
independent of the particular foliation within this subset. When the geometry 
is regular, we can place spatially invariant numerical bounds on the values of 
both the spatial and the temporal gradients of the scalar areal radius, R. 
These bounds are entirely independent of the particular gauge and of the 
magnitude of the sources. When singularities occur, we demonstrate that the 
geometry behaves in a universal way in the neighborhood of the singularity. 
These results can be exploited to develop necessary and sufficient conditions 
for the existence of both apparent horizons and singularities in the initial data 
which do not depend sensitively on the foliation. 
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I. INTRODUCTION 



In this paper we examine the constraints in general relativity when the spatial geometry 
is spherically symmetric and possesses just one asymptotically flat region This is the 

simplest gravitational scenario which exhibits local degrees of freedom. 

In [Q] we focused on the solution of the constraints when the extrinsic curvature K ab 
vanishes. Though they are simple, they nonetheless display some of the features of the 
general problem. Such solutions are, however, very special. For if the extrinsic curvature 
vanishes, the momentum constraint requires that the current density of the matter fields, 
J, must also vanish. The solutions of the Hamiltonian constraint which result therefore 
correspond to 'momentarily static' spatial geometries which do not generally occur in a 
given spacetime 0. Since, if they occur at all, they occur as isolated objects, we did not 
need to fix the foliation. Here, we extend our work to cover the general situation where 
matter flows and, as a result, the extrinsic curvature is non-vanishing. The advantage of 
having dealt separately with the momentarily static solutions is that we can focus here on 
the physical feedback on the spatial geometry introduced by extrinsic curvature. 

The introduction of extrinsic curvature complicates the analysis substantially. This oc- 
curs on two levels. The first is purely technical: the Hamiltonian constraint gets coupled to 
the momentum constraint — we have to solve a coupled system of equations. The second 
is conceptual: the constraints do not single out a unique slice through the spacetime — we 
need to specify some foliation. 

In general, the initial data is given by specifying the intrinsic and the extrinsic geome- 
try on some spacelike hypersurface which satisfy the constraints. A spherically symmetric 
geometry is completely characterized by specifying the areal radius R as a function of the 
proper radius, £. The extrinsic curvature can be expressed in a form consistent with spherical 
symmetry 

K ab = n a n b K c + (g ab - n a n b )K R , (1) 

where Kc and Kr are two spatial scalars and n a is the outward pointing unit normal to the 
two-sphere of fixed radius in the slice. 

How does one go about fixing the foliation? In principle, any foliation admitting globally 
regular solutions of the constraints is as good as any other. Ideally, therefore, we would like 
to consider a completely general slice through the spacetime; realistically, however, if the 
slicing is too general, it becomes very difficult to prove anything. At very least, the feedback 
on the spatial geometry introduced by extrinsic curvature should reflect the strength of 
the material currents flowing — the gauge certainly should not overshadow completely the 
underlying physics. 

We will focus in this paper on an 'extrinsic time' foliation of spacetime. This involves 
fixing some spatial scalar function of the extrinsic curvature. For a a spherically symmetric 
geometry, we can cast this relationship in terms of two scalars appearing in (|IJ) as follows 

F[K R ,K c ]=0, (2) 

with a possible dependence on R and £ which we have not indicated explicitly. Any gauge of 
this form should (at least implicitly) be solved to fix one of the scalars [Kc, sav ) appearing 
in (|l|), in terms of the other. 
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All previous work on the constraints in spherically symmetric relativity has focused 
exclusively on some given foliation of this type. These have been maximal slicing or the 
so called polar slicing |4j]. The latter slicing mimics the K a j, = form of the Hamiltonian 
constraint and is the foliation which provides the standard presentation of the Schwarzschild 
geometry. Unfortunately, in either case, one is at a loss to know just how sensitively the 
solution depends on the choice of gauge. How will our notions of the size and the energy 
content change in another foliation? If they change in a way we cannot quantify they are 
almost useless. To address this kind of question it is desirable to work with as large a class 
of foliations as possible. 

In [|TJ, we introduced a function a, the ratio of the two scalars defining the extrinsic 
curvature 

K c + aK R = 0. (3) 

By setting a equal to some specified function, a = a[Kn, R, £] say, Eq.@ defines an extrinsic 
time foliation. If a is a function only of R and I (in particular, if it is constant), the 
momentum constraint is exactly solvable for Kr. 

It was shown in [|TJ that each constant value of a which is greater than 0.5 provides a 
globally regular slice for appropriate sources. These values of a correspond to a spacelike 
extrinsic curvature 'vector' with respect to the superspace metric. As special cases we recover 
both the maximal slicing, with a = 2, and the polar gauge when a — > oo. 

Remarkably, one can show that even when a is not a constant the gauge continues to 
provide regular slices of spacetime so long as a < 0.5 asymptotically. All spacelike extrinsic 
curvatures in superspace provide regular foliations. The identification of potentially singular 
geometries will, however, require that the gradients of a be appropriately bounded. 

The gauges we consider, in fact, represent a very large class of extrinsic time foliations. 
Recasting (^]) in the homogeneous form, Kjc = —a[KR,R,£]Kii, ensures that when the 
material current J vanishes, = 0. In particular, flat spacetime will be foliated by 
flat spatial hypersurfaces. Indeed, when the momentum constraint is satisfied, the extrinsic 
curvature is (quasi-) linear in J, albeit in a non-local way. In this way the extrinsic curvature 
of the hypersurface responds directly to the movement of matter on it - - a physically 
reasonable criterion. 

Within this large class of extrinsic time foliations, there are universal properties exhibited 
by solutions of the constraints which are either independent of, or do not depend sensitively 
on the particular foliation. These properties divide naturally into those of globally regular 
geometries and those of singular geometries. 

In [|J, we examined these properties when the initial data was momentarily static. We 
first identified a geometrical bound on the spatial gradient of the areal radius R (the prime 
is a derivative with respect to proper radius), 

- 1 < R' < 1 , (4) 

independent of the source which was valid in all globally regular geometries. This bound 
was seen to operate at a more fundamental level than the positivity of the ADM mass. We 
then went on to investigate how the matter content of the slice can potentially force the 
appearance of either apparent horizons or singularities. We showed that when singularities 



3 



occurred, they possessed a universal form and we could place bounds on the rate of divergence 
of geometrical scalars. 

How do these results generalize? 

In globally regular geometries, the spacetime gradients of R are bounded. Firstly, when 
the weak energy condition is satisfied and a > 0.5 everywhere, the bound Eq.(^j) on the 
spatial gradient continues to hold. Secondly, and perhaps more surprising, an analogous 
bound can be placed on the extrinsic curvature when the dominant energy condition is 
satisfied. We obtain the highly non-trivial result that 

- 1 < RK R < 1 , (5) 

if a > 1. The bound ([5]) can be interpreted as a bound on R, the derivative of R with 
respect to normal proper timeBoth of these bounds are independent of the source magnitude. 
They will play a central role in the establishment of sufficient conditions for the appearance 
of apparent horizons and singularities ||. If a < 1, no such bound exists — indeed, a 
counterexample can be constructed. 

Singular geometries can occur even though both p and J are finite. The only way that 
the geometry can become singular, however, is by pinching off at some finite proper radius 
from the center. Generically, at this radius (is say), R will vanish non-analyically 

R~C(£ s -t)tt . (6) 

Remarkably, the quasi-local mass (QLM) remains finite even when the geometry is singular. 
Indeed, we show that this is always true regardless of the gauge condition. Our ability to 
identify universal behavior of this form will be crucial for the establishment of necessary 
conditions for singular geometries in a subsequent publication ||. 

Generally, the singularities of the three-geometry consistent with the constraints will be 
more severe than those which are admissable at a moment of time symmetry. If, however, 
the movement of matter is tuned so that the extrinsic curvature vanishes as the singularity is 
approached, the strength of the singularity will be determined entirely by the QLM, exactly 
as it is at a moment of time symmetry ||. We show that this tuning corresponds to an 
integrability condition on the current. If, in addition, the tuning is refined so that the QLM 
also vanishes as we approach the singularity the curvature singularity disappears and the 
spatial geometry pinches off in a regular way. This latter integrability condition involving 
the QLM is completely analogous to the integrability condition encountered at a moment of 
time symmetry. Regularity at the singularity is, of course, precisely the condition that the 
interior be a regular closed universe. If the matter fields carry conserved charges these will, 
in their turn, have integrability conditions associated with them. Viewed this way, regular 
closed universes appear to be very special universes f?|. 

The paper is organized as follows: 

We begin in Sect. 2 with a discussion of the solution of the momentum constraint. In 
Sect. 3 we provide a derivation of the bounds on R' and RKr. In Sect. 4, we derive Eq.(|). 
In Sect. 5, we discuss the integrability conditions and comment briefly on the regularity of 
Euclidean relativity. We conclude in Sect. 6 with a brief discussion and outline of future 
work. 
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II. THE CONSTRAINTS 



In this section we examine the analytical structure of the constraints when K a j, 7^ 0. 
We recall that the constraints can be written as 

K R [K R + 2K C ] - ± [2 (RR')' - R' 2 - l] = 8np (7) 

and 

K' a +£{K R -Kc) = 4nJ, (8) 

where the line element on the spatial geometry is parametrized by 

ds 2 = df + R 2 dtt 2 , (9) 

£ is the proper radial distance on the hypersurface, R is the areal radius, and we have 
expanded the extrinsic curvature according to Eq. ([!]). All derivatives are with respect to 
the proper radius of the spherical geometry, £. The spatial geometries we consider consist of 
a single asymptotically flat region with a regular center, £ = 0. We will subsequently refer 
to such geometries as regular. The appropriate boundary condition on the metric at £ = 
is then 

R(0)=0. (10) 

We recall that R'(0) = 1 if the geometry is regular at this point. We assume that both p 
and J are appropriately bounded functions of £ on some compact support. This compact 
support restriction could be easily relaxed so as to consider solutions where both p and J 
decay appropriately as one approaches infinity with little extra effort but also with little 
extra insight. 

We define the quasi-local mass m as follows 

m=!(l-J2«) + (11) 

When the constraints Eq.(|7|) and (H) and the boundary condition at the origin, Eq . fllTf ) are 
satisfied, m is determined by the sources as follows: 



m = 4vr / d£R 2 



pR' + JRK R \ . (12) 

This way, m arises as a first integral of the constraints. These equations are gauge invariant. 
In a globally regular geometry, m coincides at infinity with the ADM mass, 171^. As we 
found in || in a simpler context, the introduction of m is extremely useful and will be 
exploited repeatedly in our analysis. 

To solve the constraints classically, we need to specify some foliation. In this paper, we 
will focus on a gauge condition of the general form (^) where a is some specified function of 
the configuration variables, R, K R and £. It is possible to provide a geometrical interpretation 
for these gauges. To begin with, we know that when a = 2 this condition specifies a spacelike 
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hypersurface with maximum volume in spacetime: the trace of the extrinsic curvature K = 
K C + 2K R vanishes. It is simple to show that when a = a(£), Eq.(§) is precisely the condition 
that the modified spatial volume of a closed ball, 

V a = Aw J o d£ R a , (13) 

be a maximum. 

When a is a constant, the momentum constraint can be solved uniquely for K R in terms 
of the radial flow of matter, J, as follows 

Air r l 

K R = — — / d£R l+a J, (14) 
H R 1+a Jo v ; 

where we have exploited the regularity of the geometry at the origin to set K R (0) = 0. 

When Eqs.([|) and ( ]T4| ) are substituted into Eq.(^), we obtain a second order singular 
non- linear integro-ODE for R ||. Subject to the boundary condition, ([H]), the solution 
is uniquely determined. Not only is the extrinsic curvature completely determined by the 
material sources, so also is the spatial geometry. There are no independent gravitational 
degrees of freedom, exactly as expected. We note that in the gauge Eq.(|3|), the spatial 
geometry does not depend on the global sign of J. Of course, if we reverse the sign of J, 
the extrinsic curvature picks up a negative sign. 

When a is not constant, it is still possible to mimic the solution when a is constant. To 
do this we recast the momentum constraint in the form 

(R 1+a K R )' = A7cR 1+a J + a' ln(R/ L)R 1+a K R . (15) 

where L is any characteristic length scale. The spatial variation of a has been absorbed 
completely into the second term on the RHS. The solution is given by 

Kr= r^L di ^ R ja (^ £ ) > ( 16 ) 

where 

A(£ u £) = e^ d£2a ' HR/L) . (17) 

The constant a result is simply modulated by an exponential multiplicative correction, A. 
We note in particular that K R = when J = 0. This provides a very strong justification 
for casting the gauge in the form, (||). If a = a(R,£) alone, K R will also be linear and 
homogeneous in J. If we admit a K R dependence explicitly into a the linear correlation of 
K R and J no longer holds. K R will nonetheless be positive when J is. If, however, a is an 
even function of K R then K R will echo the parity of J: K R [— J] = —K R [J}. 

It is clear from inspection of the definition of A that spatial variations in a are anti- 
screened: remote source contributions get distorted more than nearby ones. This is poten- 
tially worrysome but, as we will see, it is not a serious obstacle. 

We will now look more closely at the analytic structure of solutions. 

Let us first focus on regular geometries. 
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III. GLOBALLY REGULAR GEOMETRIES 



We first comment on the behavior of the spatial geometry in the neighborhood of the 
origin. In fact, in the neighborhood of £ = 0, we have R ~ £, so that || 

Thus if Kr is regular at the origin then it must also vanish there. This is not, however, 
surprising. Spherical symmetry is very restrictive leaving a regular geometry no freedom to 
evolve at the origin. 

We can expand R(£) in a power series in the neighbourhood of £ = and substitute into 
Eq.® to get 

Air 

R(£)=£+ Y p(0)f + .... (19) 

A consequence of the vanishing of Kr(0) is that J will only show up at order five in this ex- 
pansion — two orders behind p. The metric at the origin clearly is not sensitively dependent 
on the current flowing there. 

The other region we need to check is outside the source. What constraints does asymp- 
totic flatness place on a? To recover an asymptotically flat spatial geometry we require 
that R(£) ~ £ to leading order. We have that K R ~ constant A (£ , £)/R 1+a . Now for 
an appropriate falloff (faster than £~ l ) on a', A(£q,£) will always saturate so that we can 
absorb it into the constant. On one hand, we note that outside the source, the integral 
identity flT2| ) implies m is a constant. However, the contribution of extrinsic curvature to m 
(Eq.([TT|)) ~ l/i? 2Qoo_1 . There represents an inconsistency if a tends asymptotically to any 
value, aoo, lower than 0.5. If > 0.5, not only Kr but also its contribution to m vanish 
asymptotically. Such a choice is simultaneously regular at the origin. 

We note that with strict inequality, > 0.5, m be dominated asymptotically by the 



first term in (|TT|) so that is encoded completely in the intrinsic geometry. In the limiting 
case, a = 0.5, the intrinsic and the extrinsic geometries share the burden. However, such a 
falloff invalidates the traditional expressions for the ADM mass. 
What can we say in general about globally regular geometries? 

We will demonstrate that they possess the remarkable property that for an appropriate 
dense subset of extrinsic time foliations both the spatial and temporal gradients of R are 
bounded numerically in a way which is entirely independent of the material sources and of 
a:. 

Suppose p satisfies the weak energy condition, p > 0. Consider any foliation satisfying 
Eq.(^|) with a > 0.5 everywhere. Then, if the geometry is regular, R' 2 < 1 everywhere. The 
proof is very simple and was given in I. We repeat it here to emphasise that the spatial 
variation of a does not enter: R' must be bounded in any regular geometry. We note that 
i?'(0) = 1 and R' — > 1 at infinity. Thus R' must possess an interior critical point. At this 
point R" = 0. In the gauge (|3|), the Hamiltonian constraint, Eq.(|7]) now reads at this point 

R' 2 = 1 - 87rpi? 2 - (2a - l)R 2 K 2 R . (20) 
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Both the second and the third terms on the RHS are negative. The result follows immedi- 
ately. 

A simple corollary of this result is the positivity of m. 

We can obtain analogous bounds on K R in the gauge Eq.(||). However, it is important 
to stress that without some control over J we should not expect Kr to be bounded. Let us 
therefore suppose that the dominant energy condition (DEC), | J\ < p is satisfied everywhere. 
Our experience in U suggests that when the DEC holds, the appropriate variables are the 
optical scalars defined by UlQ| , 

u ± = 2(#±RK R ), (21) 

which are (R times) the null expansions in the out-future and out-past directions. They are 
a useful set of variables to exploit when we are interested in identifying apparent horizons 
p|JT0|j5|]. The optical scalar which marks the presence of a future trapped surface is u + : 
u) + = at a future apparent horizon. 

It was shown in |TtJ (and again in JTJ) that when the dominant energy condition, \J\ < p, 
is satisfied 

M <k + ^\k\ 2 + A, (22) 

where k = Ma,x\RK\ and K is the trace of the extrinsic curvature. These bounds depend 
on the sources only through K. When K — (a = 2) they become numerical bounds which 
are independent of the sources, 

- 2 < uj± < 2 . (23) 

It then follows that 



RK R = (u + -u_)/4, (24) 

and R' = (uj + + u)-)/4 are bounded, \RK R \ < 1 and \R'\ < 1. 

Here, we would like to show that even when a ^ 2 in the gauge Eq. (|3|) it is still possible 
to bound u± by Eq. ([23|) . One way to do this is to bootstrap on Eq . (|22D . This way, one can 
bound |u;±| when a lies within the range 1 < a < 3. However, the bound does depend on a 
and diverges at the points, a = 1 and a = 3. However, one can do better. 

We showed in JTJ that it is possible to add and subtract the two constraints (|7|) and (g) 
to obtain simple equations for the spatial derivative of u + and 

(u ± )' = -8nR(p TJ)-^ (uj+uJ- - 4) ± lu + K c , (25) 
We now exploit the gauge condition, ([5]) and Q2"4] ) to recast these equations in the form, 

( w± )' = -8tcR{ P =p J) + — ((a - 1)uj + uj_ + 4 - aul) . (26) 

AR ^ ' 



Let us establish the inequality (23) for oo + . The argument is very similar to the one we used 
above to derive the bound R' 2 < 1. We note that uj + must be bounded in a regular geometry 
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and that uj + (0) = 2 and u + — > 2 at infinity. Thus uj + must possess an interior critial point. 
At this point uo' + = so that 

(a - 1)uj + uj_ +A-auj\ = 32irR 2 (p - J) (27) 

the right hand side of which is positive by hypothesis. Thus 

auj 2 + < (a - 1 V+w_ + 4 . (28) 

However, the quasilocal mass is positive, or equivalently |T|JT0| 

lu+lu- < 4 , (29) 

so that when a > 1, u> + satisfies Eq. (p3|) . 

In H , we pointed out that the 'binding energy' M — m of a regular spherically symetric 
system satisfying the dominant energy condition is positive when the slicing is maximal. In 
general, we have 

M - m = 4vr / d£R 2 Up + J)(2 - u + ) + (p- J) (2 + , (30) 
Jo 

which is manifestly positive whenever |u;±| < 2. Thus this result is also true for all values 
of a > 1. 

It is clear from Eq. (|28|) that Eq. (|23| ) cannot be extended to a < 1 — cj_u;+ is not bounded 
from below. To obtain a bound we need to exclude both future and past trapped surfaces 
so that u + u_ is positive. We then have for all positive a, 

K| <2Max(2/v^ Min ,l). (31) 

The results for u;_ is identical. 
This implies the bound on Kr. 

\RK r \ < Max(l/V^, l) ■ (32) 

We already have R' 2 < 1 for all a when the weak energy condition is satisfied. 

One can, in fact, easily construct a counterexample demonstrating explicitly that we 
should not expect to do better than Eq. (|3"T|) when a < 1. We do this by examining the 
values assumed by the optical scalars in the neighborhood of i = 0. We can combine 
Eq.dD and (|l|) to obtain 

^ ± ~2-^(p(0) T -^-J(0))£ 2 . (33) 

If the dominant energy condition is satisfied, then when a > 1, we have u± < 2 near the 
origin which is consistent with our result. If, however, a < 1 this is not necessarily the case. 
If J(0) exceeds (2 + a)p(0)/3 and a < 1, then uj + > 2 in the neighborhood of the origin. 
We note also that 

< (2 - f [P(0K 2 ]) 2 - (^) 2 -M0) 2 ^ < 4- (34) 



This is consistent with the inequality Eq. ( p9|) . Note also that the absolute maximum of the 
product oj + oj_ obtains at the boundary values t = and £ = oo and it is also the flat space 
value. When K = 0, this is also true of both uj + and u;_. In general, the absolute maximum 
of neither need occur at these points. 



9 



IV. SINGULAR GEOMETRIES 



So far we have assumed that the geometry is regular everywhere. A non-singular asymp- 
totically flat solution defined for all £ > will not, however, always exist for every specifi- 
cation of p and J. In this section, our task will be to understand what can go wrong and to 
identify the mechanism driving the geometry into a singularity. 

In Sect. 3 we showed that R' 2 < 1 in any globally regular geometry. Thus if R' 2 > 1 
anywhere the geometry must be singular. 

Let us suppose that R' 2 > 1 at some point. Then, when K a b satisfies Eq.([|), Eq.(^) 
implies 



RR» = i 



R' 2 



2 



+ ^-(l-2a)^-4vrp, (35) 



so that R" < and R' is decreasing there. This can only occur by R' falling through 
R' = —1. Once R' falls below this value it will continue decreasing monotonically thereafter. 
The surface with R' — — 1 in the configuration space therefore acts as a oneway membrane. 
Suppose that the areal radius is R when R' = —1. We know now that the solution must 
crash, i.e. R — > within a finite proper distance which is less than or equal to Rq from that 
point. 

Since R" < whenever R' = 1 we see that the surface R' — 1 in the configuration 
space also acts as a oneway membrane and the solution can only pass downwards through 
it. However, since at a regular center we have R' = 1 and R' starts to reduce as soon as we 
enter matter, it is clear that the region defined by R' > 1 is completely forbidden. 

We conclude that crashing through R' — —1 is the generic way the spatial geometry 
can become singular. Singularities with R' — — 1 at R — are also possible. They result, 
however, only for special finely tuned matter distributions. We will discuss them more fully 
below in Sect. 5. 

Putting the regular and singular results together, we have the following: if the geometry 
is globally regular, then — 1 < R' < 1 everywhere; if — 1 < R' < 1 everywhere, then the 
geometry is globally regular. 

One possible way that this method of constructing initial data for a spherically symmetric 
gravitational field can break down is that the slicing turns null. Since R is a four- dimensional 
scalar, nothing will go bad with it. On the other hand, £ is the spacelike proper distance 
along the slice and so d£ will become small and thus R' will become unboundedly large if 
the slice turns null. But we have shown that this cannot happen if we assume a > 0.5 and 
p > 0. Note that we do not have to assume that a = constant. Therefore, if we stay inside 
the lightcone of the super-metric, we stay outside the lightcone of the spacetime! For any 
slice satisfying a > 0.5 the only possible singularity is when R — > 0. 

Let us now examine more carefully the approach towards a singularity. In the neighbor- 
hood of the point £ = is > at which R = 0, Eq.(|l6D implies that 

R ~R^ ' ( ' 



where 



C a (£) = 4tt f Mx JR 1+a A(£, t x ) (37) 

J 
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is finite if A is. Kr will therefore be singular (for physically acceptable values of a) if the 
geometry pinches off unless the current is tuned such that 



C a (£ s ) = 0. (38) 

To examine the structure of singularities it is extremely useful to exploit the definition 
of the quasi-local mass introduced earlier. ^From a functional point of view, Eq.(|ll|) is 
identical to the energy integral in classical mechanics. To exploit this analogy, we therefore 
recast this equation as follows: 

i2* = 1 - *± + K 2 R R 2 , (39) 
K 



where m is given by Eq. (|T2|) and Kr by Eq. ([14]) . 

Let us suppose that m remains finite. Now, if C a does not vanish and a > 0.5, the most 
singular term in Eq. (|39D is the quadratic in Kr. This implies that 



R' 2 ~ R 2 K 2 R (40) 

in the neighborhood of R = 0, or R' 2 ~ C 2 /R 2a . Generically, therefore, R 2 diverges. The 
solution is 

R~ (^-) a+1 (is ~ ■ (41) 

If a > 0.5, such spatial singularities are more severe than the strong singularities discussed 
in [0 which are consistent with the Hamiltonian constraint at a moment of time symmetry. 
We will refer to the generic kind of singularity driven by extrinsic curvature as a strong J- 
type singularity. As a increases, the power law determining the strength of the singularity 
increases. Note that the limit a — > oo (the polar gauge discussed in I) is extremely singular. 
This is, however, a gauge artifact reflecting how poor the polar gauge really is. Unlike the 
strong singularities occurring when K a b = 0, at which the scalar curvature TZ remained 
finite, 1Z will generally blow up (just like K 2 R ~ l/(£g — £) 2 ). On dimensional grounds, we 
expect all curvature scalars to blow up as l/{£s ~ £) 2 a s we approach a singularity unless 
there is some constraint obstructing them from doing so. 

To show that the above analysis is self- consistent, we need to demonstrate that: for 
finite a', (i)the form factor A defined by Eq.flTTD remains finite; (ii) m remains finite. 

To this end, we note that 

f d£ a' In R/L < \a'\ Uax t dt \ In R/L\ . (42) 
Jo Jo 

The integrated logarithm is bounded. While the integrand diverges at R = 0, the integral 
is nevertheless well behaved. We note that 



o 



s 

ds In s = s In s — s . (43) 



Thus in particular, the form factor is well behaved on the flat solution, R = I. This is 
essentially all we need to check because as we have just seen R(£) ~ (£ — tj l / l+a ° at a 
singularity so that logarithm is multiplicatively identical to the flat space value. 
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We now confirm that m remains finite as we approach a strong singularity. We do this 
by demonstrating that the volume integral (12) is always finite. We note that for suitably 



bounded p and J, 

(pR 2 R\ JR 3 K R ) ~ (p, J)(i s - £)~(St) . (44) 

If a < 2, the integrand itself remains finite. In general, the integral will be finite if the 
exponent of (is — £) _1 is bounded by one. But (a, — 2)/(a + 1) < 1 for all finite values of 
a thus guaranteeing that the integrals over R 2 pR' and R?JKr converge. The only possible 
gap in this argument is the assumption that m remains finite. It is possible that m diverges 
fast enough when the singularity is approached so that m/R dominates R 2 K 2 R . This would 
require m to diverge faster than R l ~ 2a . The first term in Eq.(|l|) cannot give a divergent 
m as R 2 p obviously remains finite and / R'di = R in also bounded as we approach the 
singularity. Therefore we need only consider the JR?Kr term. This will diverge like R 2 ~ a . 
Let us assume that R ~ (is — i) 13 f° r some (3 > 0. Then, from Eq.fll2]), m will, at worst, 
diverge like m ~ (is~ £)( 2 ~ a )^ +1 . Now the requirement that the m term in Eq.(p9|) dominates 
implies (2 — a)/3 + 1 < (2a — l)(3 < 0. This in turn gives — (a + 1)(3 > 1, a negative (3\ 
Therefore it is clear that m(£s) is always finite. 

The sign of m will, however, depend on the details of the current flow. This is obvious 
from the definition Eq. ([□]). Even if R' 2 > 1, a sufficiently large value of Kr can render 
m positive. In particular, unlike the value of m assumed at strong p-singularities when 
Kab = (discussed in ||) which is always negative, the sign generally can assume either 
value. Indeed m need never even be negative in a singular geometry. Though R' decreases 
monotonically, m nonetheless can remain positive. There is no conflict with the positive 
QLM theorem. In our examination of momentarily static configurations in [[|, we found 
that m is positive everywhere except at the origin or in a neighborhood of it if and only if 
the geometry is non-singular. This is a consequence of the coincidence of the converse of 
the bounded R' 2 lemma and the converse of the positive QLM theorem when K a b = 0. In 
the general case, when K ab 7^ 0, no such coincidence occurs. 

The mass is finite at R = independent of the gauge 

We have shown above that the quasi-local mass m is finite even if the spatial slice is 
strongly singular with R — > so long as (p, J) remain finite. This result was derived on 
the assumption that the slice was chosen to satisfy the gauge condition, @ with a > 0.5. 
It turns out that the finiteness of m holds on any slice. To see this we need only consider 
Eq.(12), which is slicing independent, and Eq . (|39|) , which is effectively the definition of m. 



As we argued above, the first term in Eq.(|T2]) remains finite as we approach the singularity. 
We know that R 2 p is bounded and / d£ R' = R also is well behaved. Thus we need only 
focus on the / dl R 2 JRKr term. We know that R 2 J is bounded so we only need to control 
the RKr term. The necessary control is given by Eq. (|39|) . If R 2 K\ is the dominant term 
on the right hand side of Eq.(|39"|), we get that \RKr\ ~ —R' so the integral in Eq. (|i~2l) is 
finite as R — > 0. If m becomes large and negative so that the term —2m/R dominates, 
we get \RKr\ < —R' so again the integral converges. The only case left to consider is 
the possibility that 2m /R is positive and diverges at the same rate as R 2 K\ and some 
cancellation occurs so that R! is uncorrected with RKr. Let us assume R ~ (is — f-Y for 



12 



some p > and \RK R \ ~ (£ s - fp 7 for 7 > 0. We then get m ~ (£ s - £f- 2 ^ . However, the 
argument in Eq.fll^) goes like R 3 K R ~ (£3 — £) 2/3-7 • Hence we get m diverging at worst like 
m ~ (£3 — £) 2/3_7+1 . If this is self consistent, we require 2(3 — / y + l = (3 — 27. This implies 
(3 + 7 + 1 = 0, which makes no sense. Thus such cancellation cannot take place. 



V. INTEGRABILITY CONDITIONS 

What are the implications of the integrability condition, Eq. (|38|) ? If Eq.(p8|) is satisfied 
the strong J singularity is moderated to one which is only strong a la p. The behavior in the 
vicinity of the singularity will then be determined by the m/R term in Eq.(^) even if the 
system was originally 'driven' towards the singularity by extrinsic curvature. If, in addition, 

m{£ s ) = 4vr [ £S di \ P R 2 R' + JR 3 K R ] = , (45) 



the singularity will be a weak one with R'(£g) = — 1. We note that R" = at this point. The 
corresponding bag of gold will be a regular closed universe. These integrability conditions 
do depend on a. If a given function J satisfies Eq.fl3"8|) with one function a, generally it will 
not satisfy that condition with any other function. There is no spacetime diffeomorphism 
invariant statement of the integration. The integrability condition need not be preserved by 
the evolution. 

If J is positive (or negative) everywhere, C a (£) defined by Eq.(|37f) cannot vanish. Thus, 
if matter is collapsing or exploding everywhere, all singularities must be strong J-type 
singularities. This contrasts with the obstruction, p' < 0, discussed in p|, prohibiting the 
formation of any singularity when = 0. In general, we note that on performing an 
integration by parts on the first term, m can be rewritten 

m (£) = — p R 3 + 4vr / d£R 3 [JK R - p'\ . (46) 
3 Jo 

The first term is manifestly positive. So is the third if p' < 0. If J is positive (negative) 
everywhere then so is m in any a - gauge. However, the third term appearing on the 
RHS of Eq.(|39D may still pull the geometry into a singularity if J is sufficiently large. The 
peculiarity of momentarily static configurations with p' < discussed in 0] can clearly be 
destabilized by the motion of matter. All regular closed cosmologies simultaneously satisfy 
two integrability conditions, Eqs.([3^) and fl45|). There can be no net flow of material from 
one pole to the other. In particular, J must change sign between the poles. In addition, 
Eq.QIBD tells us that 

m(£ s ) = 4tt f S d£ R 3 [JK R - p'} = . (47) 
J 

In particular, JKr — p' must change sign between the poles. These conditions will be 
examined in the closed cosmological context in a subsequent publication [4]. 

There are no strong J singularities in the Euclidean Theory 

The singularity structure we have investigated has one important consequence for Eu- 
clidean general relativity. If the sign of the quadratic term in K R appearing in Eq. (|39|) had 
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been negative, instead of facilitating the occurrence of singularities it would have presented 
an obstacle to their occurrance. Any non- vanishing extrinsic curvature would therefore tend 
to stabilize the spatial geometry against singularity formation. We note that there is pre- 
cisely such a sign switch in the Hamiltonian constraint of Euclidean general relativity. The 
Bianchi identities there tell us that the solutions of the constraints represent all possible 
configurations the system may assume as it is evolved with respect to Euclidean time. This 
suggests that gravitational instantons will tend to be more regular than their Lorenzian 
counterparts. In fact, the most singular Euclidean geometries will occur when the geometry 
is momentarily static. In a tunneling Euclidean four-geometry, such three-geometries corre- 
spond to the initial and final hypersurfaces of the Lorentzian spacetimes between which it 
interpolates. If these hypersurfaces are themselves non-singular, i.e. do not involve Planck 
scale structures, then Planck Scale physics does not enter the semi-classical description of 
tunneling between them. This would appear to validate the application of the semi-classical 
approximation. 



VI. CONCLUSIONS 



We have identified a dense subset of extrinsic time foliations with respect to which there 
exist universal bounds on certain geometrical invariants. When the geometry is regular, 
we have described how the spacetime gradients of R in this dense subset are bounded 
numerically, independent both of the gauge and of the sources. Near a singularity, these 
gradients diverge in a way we can quantify. 

These results can be applied to address the question of identifying necessary and sufficient 
conditions for the presence of apparent horizons and singularities in the initial data [|5]|| 
extending the work of pT2| , |T3|| , |2| and [[10]. In the analysis of sufficient conditions for the 



appearance of trapped surfaces and singularities, first the moment-of-time-symmetry case 
was examined ||12| , || followed by maximal slices [0 |l(|- We find that, not only can we 



extend this work to constant a but to the large class of extrinsic time foliations described 
by Eq.(^) for variable a within the range 0.5 < a < oo. We also find that we can provide 
very powerful generalizations of the necessary conditions introduced in for moment of 
time symmetry initial data to general initial data. 

There are a number of interesting spherically symmetric problems we intend to pursue. A 
very satisfactory representation of regular closed solutions of the constraints can be given as 
closed bounded trjectories in the u_) plane. In this representation R plays a subsidiary 
role. These variables suggest a novel approach to the canonical quantization of spherically 



symmetric general relativity |TJj]. Indeed, constant a foliations can be exploited to provide 



a new description of the Schwarzschild solution [15 



The next stage is the examination of the classical evolution. Write down the Einstein 
equations with respect to the optical scalar variables. Can we cast the theory in Hamiltonian 
form? If the value of these variables in the analysis of the constraints is anything to go by, 
one has every reason to expect that they will throw light on the solution of the dynamical 
Einstein equations, both analytically and numerically. Indeed these variables have recently 
been exploited to establish global existence results [T^ . 
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